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Abstract 

Analyses of the cosmic microwave background (CMB) radiation maps made 
by the Wilkinson Microwave Anisotropy Probe (WMAP) have revealed anoma- 
lies not predicted by the standard inflationary cosmology. In particular, the 
power of the quadrupole moment of the CMB fluctuations is remarkably low, 
and the quadrupole and octopole moments are aligned mutually and with the 
geometry of the Solar system. It has been suggested in the literature that mi- 
crowave sky pollution by an unidentifled dust cloud in the vicinity of the Solar 
system may be the cause for these anomalies. In this paper, we simulate the 
thermal emission by clouds of spherical homogeneous particles of several mate- 
rials. Spectral constraints from the WMAP multi- wavelength data and earlier 
infrared observations on the hypothetical dust cloud are used to determine 
the dust cloud's physical characteristics. In order for its emissivity to demon- 
strate a flat, CMB-like wavelength dependence over the WMAP wavelengths 
(3 through 14 mm), and to be invisible in the infrared light, its particles must 
be macroscopic. Silicate spheres from several millimetres in size and carbona- 
ceous particles an order of magnitude smaller will suffice. According to our 
estimates of the abundance of such particles in the Zodiacal cloud and trans- 
neptunian belt, yielding the optical depths of the order of 10^^ for each cloud, 
the Solar-system dust can well contribute 10 fiK (within an order of magni- 
tude) in the microwaves. This is not only intriguingly close to the magnitude 
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of the anomalies (about 30 //K), but also alarmingly above the presently be- 
lieved magnitude of systematic biases of the WMAP results (below 5 /xK) and, 
to an even greater degree, of the future missions with higher sensitivities, e.g. 
PLANCK. 

This manuscript has been accepted for publication in The Astrophysical Journal. 
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1 Introduction 



The five-year rel ease (IHinshaw et al.l. 120091) of tlie Wilkinson Microwave Anisotropy 
Probe (WMAP, iBennett et al.l . l2003al ) data represents the current highlight of a 
project which provided a so far unprecedented amount of high precision cosmo logi- 
cal data. On small angular scales the measurement of the cosmic microwave back- 
ground (CMB) anisotropics led to the precise confirmation of the cold dark matter 
(ACDM) model of a nearl y spatially fia t univ erse, dominated by dark energy and 



non-baryonic dark matter (jSpergel et al.l . |2006| ) . Nevertheless at large scales a num- 



ber of unexpected, potentially damaging anomalous features in the CMB have been 
reported, indicating that either our current understanding of standard inflationary 
cosmology or the data processing, including foreground removal techniques are as 
yet inadequate. 

Among these anomalies, the expansion of the CMB sky in spherical harmonic 
functions resulte d in an octopole which is unusually planar and orien ted parallel to 



the quadrupole ( Tegmark et al. . 2003 : de Oliveira-Costa et al. . 20041 ) . Three of the 



four planes determined by the quadrupole and octopole are orthogonal to the echptic 
at 99.1% confidence level (CL), and the normals to these planes are aligned with 
the direction of the cosmological dipole and with the equi noxes, inconsis t ent Av ith 
Gaussian random, statistically isotropic skies at 99.8% CL (ISchwarz et al.l . |2004| ). 



The pres ence of non - Gauss i an features in the CMB tern perature fluctuations was 
reported by ICopi et al.l (120041 ) . iLand and Magueijd (120051 ) suggested that the pres- 
ence of preferred directions in the low-order multipoles also extends to higher multi- 
poles beyond the octopole. 

Of course any of these anomalies challenges the validity of the standard sce- 
nario of inflationary cosmology which predicts scale-free, statistically isotropic and 
Gaussian random CMB temperature fluctuations and uncorrelated multipoles. How- 
ever, although it is very unlikely that these features are just a statistical fluke, their 
cosmological orig in is still an open debate. Therefore, besides various kinds of poten- 



tial new physics (jHannestad and Mersini-Houghtonl . l2005l : iMoffatl. l2005l: iJaffe et al. 



2005 : iGordon et al.l. l2005l). a compact co smic topology (Ide Oliveira-Costa et al. 



2004 : lMota et all. 12004 : ICornish et al.l .l2004f) or modified inflation (|Piaol . l2005l : lLinde 



20041 : iHunt and Sarkai . 120041) and also conventional ef f ects of Galactic foreg round 



emission (lEriksen et al.l . l2004l : ISlosar and Seljakl . |2004J : iNaselsky et al.l . |2005| ) have 



been suggested as possible physical explanations. 

The significance of the correlations of the quadrupole and octopole with the 
echpt ic plane, howeyer, h i nts at the Solar s y stem as the origin of an unac counted 



bias ( ISchwarz et al.l . l2004l : ICopi et al.l . l2006l ) . IStarkman and Schwara (120051 ) specu- 
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lated that an unknown dark cloud of dust in the solar n eighbourhood may have con- 



tributed to the microwave sky surveyed by the WMAP. iFrischI ([20051) proposed that 



the interstellar d ust trapped naagneti cally in the heliosphere can possibly explain the 



CMB anomalies. iBabich et al.l (120071 ) discussed the possibility of discovering the dust 
of the trans-neptunian belt in the WMAP data. None of the investigations, however, 
gave fully convincing explanations to solve the problem of the CMB anomalies. 

This paper opens a series of publications in which we simulate the microwave 
thermal emission of dust clouds inside or in the vicinity of the Solar system, and test 
whether these clouds can be responsible for the unexplained correlations or can be 
discovered in the CMB experiment data. Here we focus on the spectral constraints 
on material and size distribution of the particles forming a cloud invisible in the 
infrared wavelengths and imitating closely the CMB emission in the microwaves. 
Such cloud would indeed contribute to the CMB maps without us knowing about it. 
The absolute photometry and spatial geometry of candidate clouds' contributions to 
the CMB maps are left for further scrutiny in follow-up studies. 

The paper is organised as follows. Section [2] deals with the general constraints 
on an unknown cloud's spectrum, from the infrared wavelengths to microwaves. It 
shows how a cloud of macroscopic particles can be noticeable in the microwaves and 
still avoid detection in the infrared light. Candidate clouds are proposed in the Solar 
system that can add non-negligible emission in the microwaves without being easily 
recognized in the infrared light. Estimates of the microwave temperatures of some 
known dust clouds are made, supporting the case for a detailed study. Constraints 
on the microwave spectra of foreground sources from the Internal Li near Combina- 



tion (I LC) maps derived from the WMAP multi-wavelength data (iBennett et al. 



2003bh are also formulated to facilitate the determination of plausible composition 



and size distribution of the cloud. Section [3] introduces the elements of the Mie the- 
ory of light scattering necessary for our calculations of the thermal emission from 
dust particles, and borrows bibliographic s ources on the optical properties of several 



chemical compositions from the database (IHenning et al.l . Il999l ). Thermal emission 



spectra of various sample and one natural dust clouds are generated in Sec. H] and 
tested against the constraints placed by the ILC maps. Conclusions are drawn in 
Sec. El 



2 The thermal emission due to a foreground cloud 
2.1 In the Infrared Wavelengths 

The first strong spectral constraint on the hypothetical dust cloud comes from the 
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fact that so far it has not been recognized at other wavelengths, particularly in the 
infrared light, where dust is normally very bright. This is true e.g. for the Zodi- 
acal thermal emission produced predominantly by dust with temperatures around 
300 K and peaking at ~ 10 fim, and even for the Galactic dust with temperatures 
around 10 K with a maximum radiance in the far infrared. However, these clouds 
are composed of particles much smaller than a millimetre in size, and in this section 
we demonstrate how a cloud of bigger particles can avoid detection in the infrared 
light and still shine considerably bright in the microwaves. 

Figure [1] compares intensities of the CMB radiation, its average anisotropy, and 
several dust clouds between wavelengths of 10 to 3 x 10^ fim. The emissivity of 
interplanetary dust ('IPD') is simp lified to unity below 100 /xm and above this 



wavelengt h fcf. iLeinert et al.l. Il983f ). Thi s is an approxima tion to e.g. astronomical 



silicate by Draine and Lee f 19841 ) and the Griin et al. ( 1985 ) model of interplanetary 
dust in which the b ulk of cross-section a r ea is comprised by the meteoroids from 10 
to 100 /im in size. iMaris and Buriganal (120071 ) used a similar law in their analysis 
of the microwave emission by interplanetary dust. The other clouds on the plot are 
hypothetical, with the constituent particles being macroscopic with respect to all 
wavelengths, i.e. having flat emissivities from the infrared to millimetre wavelength. 
Their geometrical optical depths r were set so as to provide an extra temperature of 
30 /iK in the microwaves, i.e. close to the magnitudes of the CMB quadrupole and 
octopole. 

The interplanetary dust at 300 K is the brightest source in the infrared light. 
The CMB radiation exceeds the emission from interplanetary dust above 300 /im 
wavelength, while the corresponding threshold wavelength for the CMB anisotropy 
is 600 /im. In the WMAP wavelength range (above 3 mm), the interplanetary dust 
is three or more orders of magnitude dimmer than the CMB anisotropy, explaining 
why it had not been taken into account as a serious bias. However, the steep decline 
of interplanetary dust in Fig. [T] with wavelength increase is due to its presumed 
emissivity oc above A = 100 /im. Clouds of macroscopic particles naturally do 
not follow this trend. 

Figure [1] allows one to put constraints on the temperature and optical depth of 
the unknown cloud that is seen in the microwaves but as yet has not been recognized 
in the infrared light. Obviously it must not be brighter than the interplanetary dust 
in the infrared wavelengths between 10 and 100 /im, where the zodiacal thermal 
emission is dominant. This criterion is met by all 'clouds' in the plot, and even 
a cloud of macro-meteoroids at 300 K, i.e . an equilibr ium dust temperature near 
Earth (see Eq. ([7]) in text below, or refer to lReachl . Il988l ). would not be immediately 
rejected. 
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Figure 1: Relative intensities of the Cosmic Microwave Background radiation, its 
anisotropy, and sample dust clouds. 'IPD' stands for interplanetary dust which is 
assumed to have a flat unit emissivity below 100 /im and an emissivity oc ab ove 
this wavelength, i.e. roughly that of astronomical silicate ( iDraine and Led . Il984l ) for 



10 to 100 /im-sized interpl anetary par t icles which are dominant in the meteoroid 
flux at 1 AU from the Sun (jGriin et al.l . 119851 ). when weighted by their cross-section 
area (cf. iLeinert et al.l . Il998l ). 'Clouds' have fiat unit emissivities from the infrared 
wavelengths to the microwaves, i.e. are composed of the macro-meteoroids big with 
respect to all wavelengths displayed. Their temperatures are determined by their 
heliocentric distances (~300 K at 1 AU, 200 K at 2 AU, 100 K at 9 AU, and 30 K at 
100 AU, see Eq. ([7]) in text). The geometrical optical depth r of 'clouds' is calculated 
so as to provide an excess temperature of 30 /xK in the microwaves. The optical depth 
of 'IPD' is 10"^. 



This criterion may look too bold: even if the cloud mimics really well the emission 
maps of smaller inte rplanetary d ust particles, its addition would double the total 
bright ness. Howev e r. iReachI (1l988l ) found that the infrared emission predictions based 
on the lGriin et al.l (1l985l ) model are half the IRAS observations. Thus a considerable 
degree of freed om persists in int erpreting the infrared observations using the size 
distribution of iGriin et al.l (119851 ). even before any exotic hypothesis is suggested, 
e.g. a local enhancement of dust around WMAP near the L2 point of Sun-Earth 
which has been hidden from the earlier sufficiently sensitive observatories like IRAS 
and COBE at solar elongations > 135°. 

Figure[T]can provide only an upper limit on the cloud's brightness. An estimate of 
the lower limit is more difficult to obtain sin ce our knowledge about b ig meteoroids 
away from the Earth orbit is rather poor. ISykes and Walked (119921 ) estimate the 
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optical depth of a cometary dust trail as 10~® (lOP/Tempel 2), and provide evidence 
that the constituent particles are cm-sized. At 250 K, the trail is 2.5 fiK bright in 
the microwaves. Similar estimates can be obtai ned for asteroid dust belts, if they 
are composed of particles with a unit emissivity (iLow et al.l . Il984j ). These structural 
elements are recognized as faint but distinct features in the smooth background 
zodiacal emission (r ~ 10"''). If there were broader structures — in space or on 
celestial maps — resembling the smaller interplanetary grain distribution, they could 
be easily confused with the Zodiacal thermal emission in the infrared wavelengths. 
In what follows, we present two plausible candidates, concentrated between Earth 
and Ju piter and in t he tra ns-neptunian belt. 

The iGriin et al.l (119851 ) model is based largely on the data acquired near Earth, 
most notably the lunar rock samples covered by plentiful micro-craters. They provide 
a vast amount of information for the model's meteoroid mass distribution inferred 
over 20 orders of magnitude (10"^*^ to 10^ g), however, the Moon spins, erasing 
memory of the direction of impacts, and it samples a very limited volume in the 
Solar-system space. There is just a little clue, therefore, as to where these mete- 
oroids came from and to what extent they are representative of the entire Solar 
system. Populations of particles have been known to exist which are not sampled 
by the Moon, with the compositions and sizes not resembling those of the Earth- 
crossing meteoroids. Examples include the above-mentioned asteroid dust bands and 
cometary trails. Meteoroid streams should exist along the orbits of the comets with 
the perihelia outside 1 AU from the Sun, just as they exist along the comet orbits 
crossing the Earth's orbit, revealed by meteor showers. Low number densities and 
fluxes of these particles have so far prevented their direct and remote registration 
outside the Earth-Moon system. Theoretical modeling is necessary to attempt to 
ful fill this observational gap. 

Hughes and McBride fll990h simulated the number density distribution of me- 



teoroids from short-period comets by distributing test particles uniformly in mean 
anomaly along each comet's orbit. Most of these particles never cross the Earth's 
orbit. Their number density peaks at 2-2.5 AU and exceeds its modest near-Earth 
level up to 5.5 AU (Fig. [2]). The near-Earth level is about five times lower than the 
maximum. If their parent comets were the only source of big meteoroids (> 100 /im 
in size) reaching Earth and Moon, t hey would constitu te about 10% of the total 
cross-section area of meteoroids in the iGriin et al.l (1l985l ) model of the fiux at 1 AU 
(see also Fig. [6]in Sec. 13.31) . As the total optical depth of the zodiacal cloud is ~ 10~^, 
and if its size distribution were spatially homogeneous, then these particles would 

contribute an optical depth of just ~ 10"*^.^ 

However, the Poynting- Robertson drag ( Wyatt and Whipple . 1950l : Leinert et al. 
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19831 : iGor'kavyi et all 119971 ) causes a drift of small meteoroids toward the Sun and 
a depletion o f these particles with helioc entric distance oc or steeper. The 
simulation by iHughes and McBridd (1l990l ) shows that big meteoroids do not follow 
this trend. In their model, assuming Keplerian motion, the large particles reside 
near the parent comet orbits. Although this assumption is only accurate if the 
planetary perturbations are ignored, by including a large ensemble of comets (135 in 



Hughes and McBridd . Il990l ) one can initially account for the planetary perturbations 



already imprinted in the cometary orbit distribution (ignoring the observational se- 
lection effects on comets). More sophisticated modeling of the sources and evolution 
of meteo roids taking the Poynting-Robertson effect and mutual collisions into ac- 
count by Ishimotol ( 2000 ) independently confirms that the number of big meteoroids 
is indeed considerably higher outside 1 AU. 

Now let us try to estimate the microwave emission from the big particles based on 
Fig. [2j The total optical depth of zodiacal cloud is based on measurements confined to 
the vicinity of Earth's orbit: the total thermal emission as well as visual light reflected 
by an interplanetary meteoroid illuminated by the Sun decreases with heliocentric 
distance R sharply: oc (i.e., proportional to the incident radiation flux). The 
number density of small meteoroids (less than ~ 100 /im in siz e) dominating the 
cross-section area of the zodiacal cloud according to the model by lGriin et al.l (119851 ) 
also decreases oc R~^ or steeper. Thus three quarters of the total thermal emission 
observed e.g. in the anti-solar direction should come from within 1 AU from the Earth. 
The emission -weighted cumulative numb ers in Fig. [2] show that for the cometary 
meteoroids of Hughes and McBridel ( 1990 ) one half of the radiation would come from 
within 1 AU outside Earth's orbit (divide the emission-weighted cumulative number 
at 10 AU by that at 2 AU): their number density initially grows with distance rather 
than decays. 

In the microwaves, the thermal emission from dust particles is in the Rayleigh- 
Jeans regime and is proportional to tem peratur e . Th e temperature decrease with 

-1/2 



distance is rather slow, oc R 



[see e.g. 



Reach . 1988, or refer to the next section 



of this paper for a derivation). Thus rn o re dis tant particles are revealed, including 
virtually a ll of the iHughes and McBride (Il990h meteoroids, whereas the small parti- 
cles of the iGriin et al.l (119851 ) model are still depleted at longer distances due to the 
Poynting-Robertson drag. 

If indeed the optical depth of the big meteoroids were ~ 10~^ for the observations 
sensitive to the emis s ion fr om dust within 1 AU from the Earth, the model by 



Hughes and McBridd ( 119901 ) predicts that all meteoroids from short-period comets 
(including those beyond 2 AU from the Sun) would have a depth of about 5 x 10^*^. A 
scaling factor of 5 stems from the total cumulative number plot in Fig. [2] which shows 
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Figure 2: Radial number densities (top) and cumulative numbers (bottom) of me- 
teoroids from short-period comets on the anti-solar line of sight from Earth (cf. 



Hughes and McBridd . Il990l ). All comet orbits were populated by an equal number 
of particles distributed uniformly in mean anomaly. The cumulative numbers use 
the number density as it is (upper curve), with a normalization to 1 at 10 AU, then 
the density weighted by the equihbrium meteoroid tem perature (mi ddle curve) pro- 
portional to the inverse square root of the distance (see iReachl . Il988l . or text below), 
and then the density weighted by the total emission (lower curve) proportional to the 
inverse square distance. An observation sensitive to the total emission (in visual or 
infrared light) is inefficient to discover meteoroids along the short-period comet or- 
bits. A temperature-sensitive observation (in the microwaves, where the meteoroids 
emit in the Ray leigh- Jeans regime with an intensity proportional to temperature) 
reveals remote dust clouds much better. 
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that only 1/5 of the Hughes and McBride f 1990l ) particles are located within 1 AU 
from the Earth. However, even some of these particles will be hidden to an Earth- 
bound microwave observer due to a slow temperature decrease. The temperature- 
weighted cumulative number restricts the scaling factor to three for the microwave 
range of wavelengths. 

Note that the base value of ~ 10^*^ for the total optical depth of big mete- 
oroids with i n 1 AU from Earth in the anti-solar direction is itself derived from the 
Griin et al. (Il985h_ model adjusted to the meteoroid flux at 1 AU from the Sun 



strictly. In the iHughes and McBridd (jl990l ) model, the number density of the big 
meteoroids is two to three times higher, on average, between 1 and 2 AU from the 
Sun, than at the Earth orbit (Fig. [2]). This allows one to raise the base value ac- 
cordingly and to come to the total temperature- weighted depths between 6x10^® 
and 9xl0~^. This is in fact comparable with the visual optical depth of 10~^ of the 
zodiacal cloud! At 150 to 300 K, these meteoroids could add 9 to 27 /iK emission 
in the microwaves, provided that they have a flat unit emissivity from the infrared 
wavelengths to microwaves, and without being immediately resolvable in infrared ra- 
diation (and very likely in visual light too, dependin g though on the particle albedo). 

The assumption of iHughes and McBridd (jl990h of an equal meteoroid number 
per comet orbit can be challenged. The production rates tend to be higher for low- 
perihelion comets and the lifetimes tend to be longer for high-aphelion comets. One 
can thus argue that if a few active comets, such as IP/Halley or 2P/Encke, produce 
more dust than all others, or some particular comet orbits allow for longer survival 
times, this in fact will only raise our estimate. We used an implicit assumption that 
the emission map of the big meteoroids is broadly as smooth as that of the zodiacal 
cloud. If a single bright comet is responsible, then the distribution is strongly biased 
towards its aphelion where the particles in Keplerian orbits stay much longer. They 
would produce a smaller and brighter spot on the sky. This is especially relevant to 
the Halley-type, long-period, and Kreutzer-group comets ("sun-grazers"). 

There is another candidate cloud never observed in the microwaves: millimetre- 
sized meteoroi ds in the trans-nept unian belt. Based on thermal emission models of 
debris disks of iKrivov et al.l (120081 ). we have estimated the microwave temperatures, 
although an uncertainty of one order of magnitude should be borne in mind. Table [1] 
lists the optical depths of the trans-neptunian belt particles of four size ranges. The 
particles which are macroscopic with respect to the microwaves (from 1 cm and above 
in size) have a total depth of 2 x 10"''. At a temperature of ~ 50 K, they would add 
~ 10 /iK if their absorption efficiency were Qabs ~ 1- When the particles of all sizes 
above 1 mm are included, the total extra temperature is as high as ~ 35 fiK. In the 
infrared light, however, the trans-neptunian belt is far too dim: taking all dust grains 
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Table 1: Geometrical optical depths of the trans-neptunian belt for an observer 
loc ated in the inner S olar system. Estimates are based o n debris disk models 
by iKrivov et al.l (120081 ) and the total belt mass of 0.02 (IFuentes and Holmanl . 
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bigger than 10 /im into consideration, we obtain a total optical depth of ~ 10^^, i.e. 
two orders of magnitude higher than that of the zodiacal cloud. However, the total 
blackbody emission at ~ 30 AU from the Sun is three orders of magnitude smaller 
than near 1 AU(see Eq. ([7]) below), hence the total thermal emission from the trans- 
neptunian belt seen from Earth in the anti-solar direction is no more than 10% of 
that of the zodiacal cloud. 

Even if all or some of these particles are not the reason for the WMAP anomaly, 
e.g. due to their emissivities below unity, they cannot be totally disregarded, espe- 
cially in the next generation CMB experiments. 

2.2 In the Microwaves 

Another spectral constraint comes from the microwave observations with WMAP 
which do not reveal any significant unknown foreground. As the CMB fluctuations 
are studied in terms of the effective temperature rather than radiance, let us first de- 
scribe the conversion details, and introduce the notation of the temperature spectrum 
of a dust cloud. 

A dust cloud of the temperature To and the column absorption area a{\) at a 
given wavelength A adds a specific energy flux BxiTi)) a{X) to the cosmic microwave 
background radiation, which is at the average temperature of T = 2.725 K, given by 
Bx(T). In the CMB studies, their sum is usually described by the excess temperature 
ATx of an imaginary blackbody emission source that would emit the same flux, i.e. 
Bx{T + ATx) = B\{T) + B\(Tj)) cr(A). When solved for ATx, this equation provides 
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one with a handy formula 



AT. 



1 + 



Ea(T)E,(Td) 



-T, 



(1) 



where E^iT) = exp{hc/kT\) 



and Boltzmann constants (cf. iFinkbeiner et al.l . Il999 



1 comes from the Planck law, h and k are the Planck 

who used an approximation 



working well only in the far infrared and microwaves). 

In the Ray leigh- Jeans regime, i.e. when A ^ hc/kT at long wavelength, the 
spectral radiance B\{T) is approximated by R\{T) = 2kT/cX^ and a simple linear 
scaling is permitted: ATx = (t(A)Td. It will be used for rough estimates of the excess 
temperatures due to dust clouds in Sec. 12.11 only. 

In general, instead of BxiT + ATx) simply ATx is given, so that the excess 
temperature can be plotted in the form of a spectrum. Temperature spectra of the 
CMB anisotropy and several known foregrounds are shown in Fig. [3l The foregrounds 
are the thermal emission from galactic dust, and free- free and synchrotron emissions 
arising, respectively, from the electron-ion scattering and acceleration of cosmic ray 
electrons in magnetic fields. It is clearly seen that the CMB anisotropy is dominant 
between 3 and 6 mm, whereas at other wavelengths either galactic dust or free-free 
and synchrotron emissions are more intensive. 

Even inside the range of its dominance from 3 to 6 mm, the CMB anisotropy 
must be cleaned of the foregrounds. Several methods have been developed to perform 
this task. One method is to obtain precise maps of foregrounds at the wavelengths 
where they dominate, refine them based on relevant physical theories , extrapolate to 



the W MAP wavelengths and subtract them from the WMAP data ( IHinshaw et al. 
20071 ). For example, the galactic dust emission was mapped using the IRAS and 



COBE observations in infrared light, cleaned from the Zodiacal dust emission by 
removing everything that is not correlated with the galactic hydrogen emission 
map, and then extrapolated to the microwaves assuming traditionally the emissivity 



law A ^ ( jSchlegel et al.l . Il998l ). This method has the caveat of not knowing every 



foreground or everything about the physics of known foregrounds. 



In contrast, the idea of the Internal Linear Combination (ILC) method (iBennett et al. 



2003bl ) is to reduce foreground and noise as far as possible by weighted linear combi- 
nations of multi-frequency data without using any special assumptions about partic- 
ular spatial structures of the foregrounds. Suppose the thermodynamic temperature 
of map i and wavelength Aj can be written as a sum of the CMB and the foreground 
temperature 

Tsky(A.) = T + AT/(A,), (2) 
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Figure 3: Excess temperatures of the CMB anisotropy and known foregrounds, rel- 
ative to the CMB average temperature (T = 2.725 K), including the thermal emis- 
sion by galactic dust, free-free emission arising from electron- ion scattering, and 
synchr otron emissiori due to the acceleration of cosmic ray electrons in magnetic 
fields ( iBennett et al.l . l2003bl ). Shown are the upper and lower limits taken from 
these authors as the three foregrounds vary across the sky. 
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the CMB maps are reconstructed by co-adding the data from the five WMAP wave- 
lengths (3.3, 5.0, 7.5, 10.0 and 13.6 mm) 



k k 

Tsky = Y.^^{T + ATfiX,)) = T + Y,w,ATfiX,), (3) 

i=l 1=1 

with the constraint that the weights Wi add to one in order to preserve only signals 
with a thermal CMB spectrum, i.e. 

E^. = l- (4) 

i=l 

Since the emission of various foregrounds like free-free, synchrotron and dust emis- 
sion shows significant spatial variation, mainly along the Galactic plane, it is clear 
that also the weights Wi cannot be held constant across the whole sky. Hence, to 
accommodate the spectral variability of the foregrounds, the entire sky has been di- 
vided into 12 separate regions, within which the weights are constant and determined 
independently by the criterion that the weights have to minimize the variance 

Var(Tsky) = Var(T) + Var (j2 w^ATf{X,)^ , (5) 

so that the infiuence of the foreground emission is suppressed down to a minimum. 
The regions have been chosen such that 10 of them cover the inner Galactic plane 
while the outer Galactic plane as well as higher Galactic latitudes are covered each 
by only one region. Since the weights are constant within each region it follows that 
also any foreground emission, especially outside the Galactic plane, should have a 
constant spectrum throughout the region. Therefore, if this ILC assumption of a 
constant spectrum turns out to be inappropriate, e.g. because of an unknown dust 
cloud, the segmentation of the sky, especially outside the Galactic plane, has to be 
improved which could therefore also lea d to a modifica t ion o f the final CMB map. 



Still, a comparison of the ILC and (IHinshaw et al.l . 120071 ) low-i multipoles does 



not reveal any major disagreement. If the microwave spectrum of the unknown 
dust cloud had significant variagations, then it wou ld have been removed , at least 



partially, by the ILC method and left intact by the iHinshaw et al.l (120071 ) method, 
leading to substantial divergence between their results. The only way the cloud can 
stay hidden for both is if it has a nearly fiat spectrum close to that of the CMB. 

This leads us to the conclusion that if the anomalous multipoles are indeed caused 
by an unaccounted dust cloud, or if there is any cloud that adds an undesirable 
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foreground, the cloud must have a relatively flat temperature spectrum between 3 
and 14 mm. We return to a quantitative assessment of how flat the temperature 
spectrum needs to be, and what it implies for the cloud's particle composition and 
sizes, in Sec. IH 



3 The thermal emission by homogeneous spheri- 
cal grains 

3.1 Excerpts from the theory of thermal emission 

Dust particles absorb, scatter and emit electromagnetic radiation. Absorption and 
scattering lead to extinction of light from sources behind the particle, with the ab- 
sorbed light heating up the particle. Also, scattering and emission add new light to 
the observer's line of sight, with the radiative energy being transferred from other 
directions and wavelengths, respectively. The absorption efficiency Qabs, scattering 
efficiency Qsca and extinction efficiency Qext = Qabs + Qsca are the ratios of the 
corresponding effective cross-section area to the geometric cross-section area of the 
particle. Knowing the radius a of a spherical particle, one can easily calculate the en- 
ergy absorbed, scattered and removed from the radiation flux at a given wavelength. 

The thermal emission of light by a particle is described in a more elaborate way 
using the thermal equilibrium equation: 

/ 7ra2Q,bs(a,A)^(A)dA = 4 / vra^g.bsK A)5a(Td) dA, (6) 
Jo Jo 

where A is the wavelength, JF is the incident radiance flux, BxiTj:,) is the blackbody 
radiance at the dust particle's temperature T^. As the left-hand side provides the 
total energy absorbed from a mono-directional incident flux, the right-hand side gives 
the total energy emitted, omni-directionally. In the Solar system, where the solar 
radiation is dominant, the incident radiance flux JF can be safely replaced by the 
solar spectrum, which is close to a blackbody radiating at 5700 K, reduced at the 
distance R from the Sun by the factor {Rq/ where Rq is the Solar radius. 

By denoting the absorption efficiency averaged over the Solar spectrum with Q©, 
and the same quantity averaged over a blackbody spectrum at temperature Td 
with Q(Td), then using the Stefan-Boltzm ann law and the Solar constant, one can 
rewrite Eq. ([6]) in a more concise form (cf. Reachl . 1988 ): 



Td = 2im[Q^/Q{T^)Y"R-"\ (7) 
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where R is measured in AU. A perfect black body with Q^bs = 1 throughout the 
spectrum has therefore a temperature of 279 K at 1 AU from the Sun, inversely 
proportional to the square root of the distance. This inverse-square-root trend is 
often closely followed by the real dust particles. 

The Mie light scattering theory allows one to calculate the efficiencies Qabs, Qsca 
and Qext as functions of particle radi us a and wavelength A once t he refractive index 
is provided for the particle material (IBohren and Huffmanl . 119831 ) . The index of re- 
fraction is a complex number m = n + ik, where the real part n is the inverse phase 
speed in the material with respect to the speed of light in vacuum, and the imaginary 
part k is the attenuation factor. The refractive index depends on wavelength. Under 
the assumption of homogeneity, it determines the propagation of electromagnetic 
waves inside the particle, while the assumption of spherical shape provides simpli- 
fications for determining the transformation of the waves at and near its boundary. 
The theory was formulated in the early 1900s in terms of the infinite series of spher- 
ical harmonic functions, its practical use was only made possible later on in the 
century by the development of computers. Standard Mie codes have been available 
to calculate the optical properties of particles based o n the refractive indices. 



A database of optical constants (IHenning et al.l . Il999l ) provides plentiful bib 



liographic references and tables of the refractive indices necessary to predict the 
light-scattering properties of dust particles composed of various materials, in a wide 
range of wavelengths. Even though the microwave range is covered very sparsely, 
one can find a number of directions to the relevant laboratory studies and remote 
observations. 



3.2 Particle composition and optical constants 



In our selection of materials, we followed partly the rationale of iReachI (119881 ) for 
materials constituting Solar system dust, largely supported by asteroid taxonomy. 
Carbonaceous particles are abundant in interstellar space and have been directly 
observed with the mass spectrometers on board Giotto and V ega- 2 during the ren- 
dezvous with comet IP/Halley in 1986 (jjessberger et al.l . 119881 ). Carbonaceous mate- 
rial covers the surfaces of C-type asteroids which is the dominant type in the asteroid 
belt, with 75% of known asteroids falling in this category. Silicate material largely 
constitutes the lunar and terrestrial rocks and was also revealed in the mass spectra 
of the IP/Halley dust. It is abundant in circumstellar debris disks, all o wing for re- 



mote studies of optical properties of silicate particles (lOssenkopf et al.l . Il992l ). 17% 
of known asteroids are of S-type, i.e. silicaceous according to their surface spectra. 
Even though the cosmic abundance of Fe is low, iron is the major component of some 
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meteorites, magnetite (Fe304) is present in many of them. L-type asteroids, which 
constitute 7% of known asteroid s show metallic surfaces. In more distant parts of 
the Solar system not reviewed by iReachI (1l988l ) , ice grains are an impo rtant popula- 



tion. Ice has been reported to cover the trans-neptunian objects (see iBrown et al. 



19991 : iJewitt and Luul . l2004l ). The icy moons and rings of Jupiter and Saturn further 



support the case for including ice in our review. Therefore, we will focus ourselves 
on the optical prope r ties o f carbonaceous, silicate, iron and icy particles. 



Ossenkopf et al.l (119921 ) inferred the light scattering properties of the silicate 



grains from the opacities of circumstellar dust disks at the wavelengths of up to 
10 mm. Between 1 and 10 mm, the real part of the refractive index stays near 
n ^ 3, while the imaginary part drops from k ~ 0.025 at the lower boundary to 
k ^ 0.0025 at the upper boundary of the range, inversely proportional to wavelength 
(Fig. |l^,b). This is reported both for warm oxygen-deficient circumstellar silicates 
and cool oxygen-rich interstellar silicates cons idered in the paper. 

The optical constants for silicate particle s (ILaor and Draind. Il993l . Fig. llK,b) dif- 



fer only slightly from the values given by lOssenkopf et al.l (119921 ) . with a similar 
inverse phase speed n = 3.4 and an attenuation factor decreasing with wavelength 
at the same rate yet f rom a higher k ~ 0.05 at 1 mm. Note, however, that the 
Laor and Draind (119931 ) constants are not provided for A > 1 mm and had to be ex- 



trapolated, assuming that the steep downward tr end is continu e d in t he microwaves. 

The optical co nstants for olivin e prov ided by iPoUack et al.l (119941 ). however, dis- 
agree with those of lOssenkopf et al.l (Il992l) and a simple extrap olation of lLaor and Draine 
(I1993I ) beyond a wavelength of 1 mm. lOssenkopf et al.l (1l992l ) do not quote any ob- 
servations beyond 1 mm. It is reasonable to assume that these authors extrapolated 
the trend seen in the far i nfrared wavelength into the microwaves. In contrast, 
the const ants for olivine bv iPollack et al.l (Il994j ) are based on laboratory measure- 
ments by ICampbell and Ulrichd (119691 ) at 8.57 mm and 66.7 mm who show that 
silicates do not gain t r anspa rency in the microwaves as fast as the extrapolations 



suggest. iBoudet et al.l (120051 ) may have found an explanation for the disagreement. 



They studied the temperature dependence of the absorption by amorphous silicate 
grains between 10 and 300 K and found that with increasing temperatures the ab- 
sorption efficiency grows considerably (by an order of magnitude) already in sub- 
millimetre wavelengths. Therefore , the theoretical considerations behind the extrap- 
olation bv lOssenkopf et al.l (1l992l ) may still hold true for cold dust, whereas the 
Campbell and Ulrichsl (1969) measurements in the warmth of a ground-based labo- 
ratory turn out to be more relevant to the near-Earth dust environme nt. Since ou r 



main application of the Mie theory will be for a higher temperature, the lPoUack et al. 
(j 19941 ) model is adopted as most reliable beyond 1 mm. 
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Figure 4: Optical constants fo r several che n iical co mpositions of s p herica l 
homogeneous particles found in iPoUack et al.l (Il994l ): lOssenkopf et al.l (Il992l ): 
Laor and Draind (Il993l ): IZubko et al.l (|l996l ). Top to bottom: silicate, carbonaceous, 
and iron particles. Thick curves are the data provided by the above authors, thin 
curves are our extrapolations. The solid curve in each plot emphasizes the vari- 
ant that is chosen for further calculations of the absorption efficiencies and thermal 



emission. 
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In c ontrast to sili c ates, the amorphous carbonaceous grains studied in the labora- 
tory by IZubko et al.l (1l996l ) show much higher real parts of the refractive index n of 
up to ten and more, as well as imaginary parts k of several at the wavelengths near 
1 mm (Fig. Ilt;,d). The imaginary part of the refractive index does not fall but grows 
with increasing wavelength for the 'BE' sample, i.e. amorphous carbon grains pro- 
du ced in burning benz ene in air under normal conditions. Two other species studied 
by lZubko et al.l (119961 ). those produced by arc discharge between amorphous carbon 
electrodes in different atmospheres ('AGAR' and 'ACH2'), are characterised by the 
imaginary parts of the refractive indices reaching maxima between 100 and 1000 ^m. 
and then turning down. Nevertheless, as we have checked, they show qualitatively 
similar dependencies of the absorption efficiency on wavelength. 

Additionally, the optic al constants for graphit e particles up to 1 mm wavelength 
are provided in the paper (jLaor and Draind . Il993l . Fig. |lt,d). The graphite particles 
behave similarly to the above-listed species of amorphous carbon when the electric 
field vector is perpend icular to the plane of graphite cleaves (their || c' explained in 



Draine and Led . 1 19841 ). while the particles with the cleaves being parallel to the field 
vector are significantly different {^E ± c'). The absorption efficiency in the latter case 
is very low, however, making the thermal emission by a cloud of randomly oriented 
graphite particles well described by the first case. 



We adopt for further calculations the optical constants of the 'BE' sample (jZubko et al. 



19961 ) as a reasonably good representation for the absorption efficiency of the above- 



list ed carbonaceous sp ecies. 



Zubko et al.l ( 1l996l ) have limited their study to wavelengths up to 2 mm, while 



our application of the Mie theory to the WMAP data requires the optical constants 
for up to ~ 10 mm. We have not found any constants fo r carbonaceous part icles 
at A > 2 mm in the literature. The optical constants of I Zubko et al.l (119961 ) are 
therefore extrapolated as shown in Fig. Ill:;,d (thin lines). 

Iron particle s (Fig. |l^,f ) have a rather plain dependence of the optical constants 



on wavelength (jPoUack et al.l . Il994l ). As shown in that paper, qualitatively similar 



dependencies are also exhibited by iron combinations with some other elements, 
such as iron sulfide FeS, a circumstance that allows the results of the metallic iron 
model application to be expanded onto a correspondingly broader range of particle 
composition. 

Water ice is reported to have n ~ 1.8 and k decaying fr om 10~^ at the wavelength 
of 1 mm to 10-3 at 10 mm, and at a temperature of — 1°C (IWarrenl . 1 19841 ). For lower 
temperatures, k is lower: at — 60°C, the values range from 3 x 10"^ at 1 mm to 
3 X 10""^ at 10 mm. Solid ice can survive for a long time only far from the Sun, where 
the temperatures are low (see Eq. [7]), e.g. ~ 50 K = — 223°C in the trans-neptunian 
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belt. The attenuation factors for cosmic ice can therefore be even lower than 10~^. 
Then the grains absorb or emit virtually no microwave radiation, unless they are 
much bigger than th e wavelength, i.e. centimetres in size. This is in agreement with 
Pollack et al.l (119941 ). who show that the imaginary part decays from 3-10~^ to 3-10~^, 
while the real part of the refractive index is in agreement with the value found by 
WarrenI (119841 ). It is noteworthy that CO2 ice has a similar real part n ^ 1.4, whereas 
the i maginary part of the refractive index almost vanishes (m < 10~^) (jWarrenl . 



19861 ). While the attenuation factor can be raised by impurities, it is generally fair 
to say that cosmic ice absorbs the microwave radiation negligibly with respect to 
carbons and even silicates. Therefore we discuss neither water nor CO2 ice in the 
remainder of the paper. 



3.3 Calculation of absorption efficiencies 

The absorption efficiencies for the particle compositions considered above, except for 
ice, are shown in Fig. [51 The left panels show single-particle efficiencies, while in the 
right panels the efficiencies are averaged over different size distributions that can be 
expected in dust clouds. 

Single-particle plots show a steep decrease of the absorption efficiencies of small 
grains with wavelength. Because the attenuation factor sharply weakens with wave- 
length, the absorption efficiency is proportional to Qabs for silicate particles 
of 1 Atm radius b etwee n 100 and 600 /xm wavelengths. If the optical constants by 
Ossenkopf et al.l (119921 ) were adopted, this trend wo u ld con tinue to the microwaves. 
However, the flattening attenuation of iPoUack et al.l (119941 ) changes the trend some- 
what above 600 fim. As the particle size grows, the absorption efficiency flattens 
too, with a roughly constant value for wavelengths shorter than the particle size. 

The carbonaceous particles show an analogous qualitative behavior, however, 
due to much stronger attenuation, their absorption efficiencies are flat for A less than 
~ 10a. This is clearly seen for A < 2 mm where the actual measurements have been 
used by lZubko et al.l (119961 ) to determine the optical constants. 

Interestingly, the absorption efficiencies of iron particles starting from ~ 10 /im 
in size are much less size-dependent over the wavelength range under consideration. 
They also do not completely flatten up to the biggest size considered, i.e. ~ 10 mm, 
but instead all sizes above 10 /zm share a weak dependence Qabs oc A~^/^. Even 
though this decay is extremely weak in comparison with that of like-sized silicate 
and carbonaceous particles, it leads to the absorption efficiency vanishing to almost 
10^"^ at wavelength of several millimetres, which makes iron a very unlikely candidate 
for the elusive cloud material. For the sake of brevity it will not be considered during 
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Figure 5: Absorption efficiencies of spherical homogeneous particles for the same 
chemical compositions as in Fig. HI From top to bottom: silicate, carbonaceous and 
iron particles. Left column: absorption efficiencies for single particles of different 
sizes, right column: for size distribution s n(a) da = a~'^da and the interplanetary 
meteoroid size distribution of iGriin et al.l (119851 ) . 
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subsequent deliberations. 

The size distributions assumed in our calculations include power laws n{a)da 



a '^d a with several slopes 7, and the size distribution of interplanetary meteoroids (iGriin et al. 



19851 ) which was derived from the crater size distribution on lunar rock samples 
and remains the basis for the modern met eoroid environment models (jPivind . Il993 : 
Staubach et all 119971 : iDikarev et al.l . l2005l ) . The absorption efficiencies of single par- 
ticles were averaged over the cross-section area using 

/Qabs(Q, A)n(a)aMa 
/ n{a)a'^da ' 

where the integration limits were separated sufficiently to bracket all the dust that 
can emit at a given wavelength (ti ny submicro n-sized grains are inefficient already 
in infrared radiation, so e.g. for the Griin et al. (|l985i ) distribution of interplanetary 
meteoroids in Fig. [6] a moderate growth of cross-section area in the grains with 
s < 100 nanometres is invisible; huge meter-sized boulders are too rare). 

The slopes 7 wer e picke d around the Mathis-Rumpl-Nordsieck (MRN) distri- 
bution (IMathis et al.l . 119771 ) found for the interstellar dust and which are typical 
for many other dust clouds as well, i.e. 7 = 3.5. This slope was also derived 
for t he steady-state size distribution of a cloud of colliding and disrupting parti- 
cles ( iDohnanyil . Il969l ). The slope can be modified by particle dynamics. 

Figure M shows the size distributions used in our calculations as well as one of 
their momenta, the cross-section area distribution. Shallow slopes 7 emphasize bigger 
particles in the cloud, and make the average absorp tion efficien c y acco rdingly flatter. 
The interplanetary meteoroid size distribution by iGriin et al.l (119851 ) yields a steep 
absorption efficiency since th e bulk of its cross-s ection is comprised by particles from 
10 to 100 /im in size. The I Griin et al.l (119851 ) distribution has different slopes in 
various size ranges, where different dynamics determine distinct particle lifetimes: a 
slope close to 3.5 for a < 1 fim, near 2 between 1 and 100 /xm, and about 5 above 
100 yum. 



3.4 Nonsphericity and inhomogeneity 

The assumptions of sphericity and homogeneity made above for the sake of sim- 
plicity can of course be debated. Aggregates can indeed emit differently at long 
wavelengths, especially if they include conducting materials (carbon, iron). Interest- 
ingly, by taking the porosity into account one can come to flat absorption efficiencies 
for even smaller particles. Carbonaceous ballistic particle-cluster aggregates (PCAs) 
usually give flatter efficiencies than Mie spheres. The more realistic the method. 
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Figure 6: Size distributions n(a) da = a~'^da and the interplanetary meteoroid size 
distribution of iGriin et al.l (119851 ) used in producing the average absorption efficien- 
cies in Fig. |5] (top). The corresponding cross-section area distributions (bottom). 
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the more pronounced is the effect: e.g. modified spectral function (MSF) predicts 
flatter efficiencies than discrete multipole method (DMM), and D MM flatter than 



effective-medium theories (EMT, see e.g. Stognienko et al. . 19951 . for more details 



and references). Therefore, we have good reasons to trust that our Mie calculations 
provide results rather conservative with regard to flatness, which is important for 
our conclusions. 



4 The temperature spectra of dust clouds 

The column absorption area of a dust cloud is simply the product of its column 
cross-section area and the absorption efficiency of constituting particles. Assuming 
the column cross-section area to equal 10"'', i.e. close to that of the Zodiacal cloud 
near Earth, and taking Qabs from Sec.[3l we calculated the corresponding temperature 
spectra (see Fig. [7]). The dust particle temperature was set to 300 K (cf. Eq. ([7])). 

The peak of the blackbody emission at that temperature is near 10 fim, while 
the CMB maximum is near 1 mm. This leads to very high excess temperatures due 
to dust over the CMB at the short wavelengths. The temperature spectra of small, 
micrometre-sized particles drop sharply in the far infrared and microwaves, however. 
In contrast, bigger particles show a pronounced flattening of temperature spectra 
if their size is similar to or bigger than the wavelength. The thermal emission by 
dust and the CMB is on the Rayleigh- Jeans side of their Planck functions above 
1 mm wavelength, where they are approaching proportionality to temperature, and 
the temperature spectrum of dust is determined mainly by the absorption efficiency, 
which is flat for macroscopic particles. Silicate particles above several millimetres in 
radius exhibit flat temperature spectra, while the emission from carbonaceous grains 
larger than several hundreds of micrometres in size are already unaffected and thus 
not excluded by the ILC procedure. 

It is important to note that the temperature of the dust cloud needs not neces- 
sarily be low to produce a flat temperature spectrum over the WMAP wavelengths. 
Cooling the dust particles, e.g. by moving them away to the trans-neptunian belt 
(to 50 K) or even further to the Oort cloud, will by itself not make the temperature 
spectrum flat in the microwaves. Even warm (300 K) particles can give a rather flat 
spectrum if they are big. Simultaneously, due to their broadly flat emissivity, they 
do not reveal themselves at shorter wavelengths, contrary to small micrometre-sized 
dust. Therefore the big particles may easily be outshone by more abundant small 
particles in the i nfrared, yet d ominate the microwave emission. 



A proposal by iFrischI (120051 ) that the thermal emission of interstellar dust trapped 



in the heliosphere is an explanation of the CMB anomalies can be rejected since the 
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Figure 7: Temperature spectra of dust clouds composed of silicate (top) and car- 
bonaceous (bottom) particles. Monosize clouds are displayed in the left column, size 
distributions as in Fig. [5] are in the right column. A dust temperature of 300 K 
and column cross-section area of 10~^ were assumed, similar to those of the Zodiacal 
cloud as seen from the Earth. 



grains are too small, sub-millimetre in size. No matter how far they are from the 
Sun, due to a steep decrease of their absorption efficiencies with wavelength they 
would be rejected by the ILC map construction procedure, and visible in the infrared 
wavelengths. 

Power-law size distributions of dust are characterised by flat temperature spec- 
tra only if their slope 7 is weak, so that the big particles constitute the bulk of the 



cross- section area, according to our findings about mono-size clouds. The lGriin et al. 



(119851 ) size distribution of interplanetary meteoroids gives a temperature spectrum 
that is quite steep for the ILC if the meteoroids are silicate, shallower if the particles 
are carbonaceous. Figure [8] zooms into the WMAP wavelength range. More particle 
sizes are plotted in the mono-size cloud ca se to facilitate mo re accurate assessment 
of the temperature spectra. Note that the iGriin et al.l (119851 ) model's size distribu- 
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Figure 8: Zoom of Fig. [7] into the WMAP wavelength range. 



tion, when normahsed to the zodiacal cloud visual optical depth (10"''), is brighter 
than 1 /iK up to 6 mm wavelength under the hypothesis that all meteoroids are 
carbonaceous (bottom-right plot)! At the shortest wavelength of WMAP (3 mm) it 
can be as bright as several fiK. If indeed there are five times more big meteoroids in 
the interplanetary space outside the Earth's orbit (Sec. 12.11 and Fig. [2]) than in the 
flux at 1 AU, one can get up to ~ lOyuK thermal emission in the microwaves. Of 
course, the uncertainties of the chemical compostion and size distribution impact the 
accuracy of such estimates. Similarly, the trans-neptunian belt estimated to have a 
geometrical optical depth of 2 x 10~^ in particles bigger than 1 cm, would provide 
~ 15 /iK emission if it were fully carbonaceous (Qabs ~ 0.5 in Fig. [5l cf. Sec. 12.11) . 
This is not the there are many ice particles in the belt which are not visible 

in the microwaves, however, the number is intriguingly high. 

The "color ratio" of mono-size clouds is plotted in Fig. M (left) as a function of 
size. It is the ratio of the excess temperatures of the clouds in the 13.6 to that in 
the 3.3 mm waveband, AT13.6/AT3.3. The color ratio of the carbonaceous particles 



of IZubko et al.l (119961 ) rises to close to unity for a particle size between 1 and 2 mm. 
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i.e. where reliable optical constants are still available. The color ratio for olivine 
remains low up to particle size of several mm. Obviously, in order to pass the ILC 
filter, one needs bigger silicate particles than carbonaceous o nes. Interest ii igly, the 
color ratio of all power laws and interplanetary meteoroids of iGriin et al.l (Il985l ) is 
too low, typically below 0.5. 

We have also calculated how much emission from dust can pass through the 
ILC filters, assuming that the cloud is sufficiently "transparent" to leave the filter 
weights unaffected (i.e., determined by other foregrounds with higher brightness 
and/or steeper spectra). If the cloud were not "transparent", the weights would be 
affected and a fraction of the cloud ern i ssion would be removed, leading to a difference 
between the ILC and iHinshaw et al.l (120071 ) foreground-cleaned maps, which is too 
low (Sec. 12.21) . Clouds with higher passthrough ratios are therefore more plausible, 
even though not proven, candidates for sources of large unaccounted contaminations 
in the WMAP data. Note that clouds with low passthrough ratios would affect the 
ILC weights in a way that the ratios are further lowered. 

Figure [9] (right) displays the passthrough ratio for the ILC weights derived for 
region 1 encompassing the higher galactic latitudes and most of the ecliptic plane. 
This is the largest region on the sky where one could be most eager to search for an 
unknown dust cloud since the galactic emission is minimal inside it. It is interesting 
to note that the emission from the clouds compose d of olivine pa r ticles is reduced by 
about 30% only if the particles are microscopic. lEriksen et al.l (120041 ) had already 
emphasized that the ILC procedure is not very efficient in removing the emission 
from dust, although they found it more efficient than in our test: approximately 
half of the VF-band (3.3 mm) simulated emission from dust was in their "cleaned" 
maps. They used the emissivity which is steeper than that of olivine. In our 
test, carbonaceous material allows us to reduce the emission of clouds of microscopic 
particles by half. 

The ILC passthrough test confirms the results of the color ratio test: olivine 
particles of about 10 mm in size and carbonaceous particles of roughly 1 mm in size 
reach a nearly 100% "transparency" for the bias removal procedure. However, one 
should bear in mind that monosize clouds are very rare in nature. A broad size 
distribution should allow clouds composed of smaller particles to become "transpar- 
ent" too, as the sharp peaks and drops of the ILC passthrough at big sizes would 
cancel out after an averaging. The iGriin et al.l (1l985l ) size distribution of interplane- 
tary meteoroids has ILC passthroughs of 64% and 58% for carbonaceous and olivine 
particles, respectively. If, however, the big meteoroids are under represented in the 
meteoroid flux at 1 AU, the percentages may be higher. Ignoring all meteoroids less 
than 100 /xm in size, one can come to passthroughs of 95% (carbonaceous) and 55% 
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Figure 9: "Color ratio" (left) of olivine and carbonaceous particles as a function of 
their size. The color is defined as Tia e/Ts.a, i.e. the ratio of temperatures of the 
cloud in 13.6 and 3.3 mm wavebands of the WMAP radiometers. ILC passthrough 
ratio (right), per cent, for monosize clouds of olivine and carbonaceous particles in 
region 1 (higher galactic latitudes). The ratio is a fracti on of the cloud's emis sion 



that passes through the ILC bias correction procedure (IBennett et al.l . l2003bl ). to 



the total emission. The ratio can be above 100% as well as negative due to the 
peculiarities of the cloud spectrum. 



(olivine). These numbers are directly attributable to the cloud of meteoroids from 
short-period comets discussed in Sec. 12. 1[ Note that for real non-spherical inhomo- 
geneous particles, the percentages will likely be even closer to 100% (Sec. 13. 4p . 



5 Conclusion 

We have investigated the microwave thermal emission by dust in the Solar system. 
We applied the Mie theory of light scattering by spherical homogeneous particles 
in order to characterise the thermal emission spectra of silicate, carbonaceous, iron 
and icy particles. Our study is partly motivated by the WMAP observations of the 
CMB fluctuations that revealed large-scale structures aligned with the Solar-system 
geometry which are difficult to explain by the standard inflationary cosmology. One 
possibility to produce such fluctuations is by a dust cloud inside or in the vicinity 
of the Solar system. Another motivation for this study is to assess the feasibility of 
detection of dust in and near the Solar system in the microwaves. 

We used the WMAP multi-wavelength observations and infrared surveys to con- 
strain the physical properties of particles constituting the hypothetical cloud and 
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to estimate the microwave emission by solar-system dust. We have found that only 
macroscopic, mm-sized silicate or carbonaceous grains could produce thermal emis- 
sion with a spectrum compatible with that of the CMB fluctuations. Smaller dust 
grains, as well as iron particles, emit with a spectrum that can easily be distinguished 
from the CMB. The small particles would also be so bright in the infrared light that 
they would have been seen by the relevant telescopes, whereas the big particles have a 
fiat emissivity throughout the spectrum, so that the abundant small grains outshine 
them at infrared wavelengths. 

In order to attain the fiat emissivity at the WMAP wavelengths, silicate particles 
must be several mm in size at least, whereas carbonaceous particles can be an order of 
magnitude smaller. This mak;es the carbonaceous particles the most likely candidate, 
as small grains are typica lly more likely in dust cloud s than big ones. Meteoroids 
from short-period comets ( Hughes and McBride . 1990l ) are a plausible candidate for 
such cloud. The dust needs not necessarily be cold nor remote. When the cloud is 
composed of particles with a broad size distribution, even smaller dust grains can 
contribute to an overall fiat spectru m without revealing themselves. For example, 
all particles of the iGriin et al.l (119851 ) meteoroid fiux model with sizes above 100 /im 
have an integral spectrum sufficiently fiat to pass the ILC procedure used to clean the 
WMAP observations from biases, if they are composed of carbonaceous material. The 
trans-neptunian belt is another plausible candidate. According to our preliminary 
estimates, each candidat e cloud can emit roughl y ~ 10 fiK in the microwaves. 



We found studies (see lStognienko et al.l . ll995l . for references) showing that the Mie 



theory provides in general more pessimistic estimates of compatibility between the 
dust and CMB spectra than advanced light scattering theories taking nonsphericity 
and inhomogeneity into account. It means that the real particles can be even smaller 
than those quoted above in this conclusion. 

There is a lack of measurements of the optical constants of carbonaceous particles 
above a wavelength of 2 mm. Given that this material provides the best chances so far 
to explain the WMAP anomaly, laboratory measurements of carbonaceous particles 
would be very helpful for further studies of the problem. 

In subsequent papers of this series, using the thermal emission model introduced 
here, we will test the absolute photometry and projected geometry of various known 
and hypothetical dust clouds in and near the Solar system, against the maps of CMB 
fluctuations made by WMAP. 
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